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SUMMARY 

The    first   step   in   the    development   of   a  hybrid   com- 
puter  simulation  of   the    data-based,    5    DOF,    3K-SES    program 
that  will  operate    in  real-time   has   been   completed.      A 
3   DOF    (Flat)    turn  model    for   an    approximate    and   simplified 
XR-3    craft  was    used  to   provide    the    equations    of   motion   in 
order   to  become    familiar  with    1)       the    turning   character- 
istics  of   the    CABSES,    2)    the    requirements    of   the    (NPS)    NAV 
PG  School  hybrid   computer  with   regard   to   signal    conversion, 
timing,    accuracy    and   loop    delay    due    to  both   signal    con- 
version  and    (right-hand-side)    computation   and    3)    the 
approximate    size    of  program  that    can   be    computed   on   the 
existing  hardware   of   the  NPS    facility. 

In  parallel   with   this    investigation,    the    5    DOF   data- 
based   digital   program   for   the    3K-SES    has    been    converted 
from   CDC   to    IBM  FORTRAN    and   is   now   operational    and  will 
serve    as    the    standard    for  evaluation   of   the    final   hybrid 
simulation   for   the    3K-SES. 

Results    of   the   hybrid   computer   simulation   of   the    3    DOF 
turn  model   of   the    XR-3   indicate    that   the   installation   of   the 
5    DOF  program  of    the    3K-SES    should  proceed    (in   the    next 
phase)     as    a   straight   forward   application   of   the    right-hand- 
side    computations   within      the    digital   computer   and    addi- 
tional  degrees   of    freedom  on  the  analog  computer.      The 
successful   operation   of   the   hybrid   computer   simulation  will 
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depend  upon   ability   of  the    10   year-old  hardware—in 
particular  the    digital    computer — to    compute    the   many 
computational   operations    required   for   the    right-hand- 
side   of   the   equations    of  motion   in  the   time    frame   necessary 
to   approximate   real-time   manual   control. 
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INTRODUCTION 


Reference  1  presents  the  Rohr  Marine,  Inc.  analysis 
of  the  stability  and  maneuverability  of  the  3K-SES  under 
direct  manual  control.  The  recommendations  contained  in 
Reference  1  included  a  suggestion  to  continue  the  investi- 
gation of  manual  control  using  a  more  realistic  model  of 
the    response    of  the   human   operator   to   emergency    conditions. 

If   a   real   time    simulation    can  be    developed   it    can 
provide    a  means    of   inserting   an   actual  human   operator   in 
the    loop,    as   well    as    any   selected  model   of   the   operator. 
A  hybrid   realization   of   such    a   real    time    simulator   appears 
to   be    feasible    and  would  provide    a   convenient   hands-on 
facility    for   a   variety   of  human   operator   studies,    plus 
many   other   studies    in    control    and   guidance. 

This    report  presents   the    results   of   the    first  phase 
of    a   feasibility   study,    in  which    it   is   planned   to   develop 
a   real    time    simulation   using   the   NPS    Hybrid   Computer 
which    consists    of   a   CI-5000    analog    unit    and   an   XDS-9  300 
digital   unit.      The    initial    studies    reported  here    are 
concerned  with    development   of    3    DOF  models    for   the   XR-3 
and   the    3K,    as    a  preliminary   to    the    development   of    a   5    DOF 
Hybrid  model. 


II.   FAMILIARIZATION  AND  PLANNING 


IIA.   THE  3  DEGREE  OF  FREEDOM  MODEL 

1.   INTRODUCTION 

The  development  of  any  computer  model  requires  famili- 
arity with  both  the  system  equations  and  the  computer  to  be 
used.   In  addition,  one  needs  quantitative  data  to  insert 
the  model  into  the  computer,  and  one  also  needs  typical 
performance  data  such  as  experimental  test  results  in  order 
to  validate  and  verify  the  model.   To  develop  a  hybrid  model 
of  the  3K-SES,  we  chose  to  start  by  modelling  the  XR-3, 
because  we  have  the  required  quantitative  data,  we  have  a 
good  6  DOF  digital  program,  and  we  have  good  experimental 
data,  plus  availablility  of  the  XR-3  craft  for  test  as  needed. 
We  also  decided  to  start  with  a  3  DOF  model,  for  several 
reasons : 

a)  It  is  the  simplest  model  which  will  permit  turning 
tests  that  are  verifiable. 

b)  The  complexity  of  the  initial  model  will  thus  be 
minimal,  but  we  can  readily  add  two  more  degrees 
of  freedom  with  confidence  in  the  resulting  5  DOF 
model,  because  the  similarity  of  the  differential 
equations  will  permit  ready  application  of  the 
techniques  we  develop  in  perfecting  the  3  DOF  model. 

In  addition  to  this  decision  to  start  with  a  3  DOF  model, 
we  chose  to  simplify  the  model  by  simplifying  the  right-hand- 
side  of  the  equations.   Since  the  terms  on  the  right-hand- 
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side  are  arithmetic  calculations  to  be  done  in  the  digital 
unit,  and  because  the  many  terms  are  additive ,  truncation 
of  the  right -hand- side,  with  simplification  of  the  terms 
retained,  should  have  little  effect  on  results  except  some 
decrease  in  quantitative  accuracy.   This  permits  us  to  work 
with  a  less  complex  model  in  the  initial  stages,  and  addition 
of  terms  at  a  later  point  should  be  a  simple  matter.   It 
will  also  have  the  advantage  of  permitting  us  to  assess  the 
effects  of  adding  terms  on  the  computational  delays  and  on 
the  limits  of  capability  of  our  hybrid  installation. 

2.   EQUATIONS  OF  MOTION  OF  THE  3  DOF  MODEL 

The  three  degrees  of  freedom  to  be  used  in  this  devel- 
opment are:   SURGE,  SWAY,  YAW.   As  a  result  the  following 
assumptions  and  approximations  are  needed: 

1.  Roll  is  zero 

2.  Pitch  and  draft  remain  constant 

3.  The  origin  of  the  ship's  coordinate  system  is 
located  at  the  center  of  gravity. 

4 .  At  t  =  0  the  ship  coordinate  system  and  the  fixed 
(geographical)  coordinate  system  have  the  same 
origin  and  the  same  coordinate  axes. 

5.  At  t  =  0  the  ship  is  heading  in  the  positive 

X    direction, 
nav 
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Figure  1.  DEFINITION  OF  COORDINATE  SYSTEMS. 


*NAV 


ORIGIN 
OF 
FIXED 
(GEOGRAPHIC) 
COORDINATES 


ORIGIN   OF 
SHIP'S   COORDINATES. 


SHIPs  PATH 


K*> 


CLOCKWISE 
ROTATION  IS  PQi 
ANGLE   DIRECTIO 


NOTE  ■•     RUDDER  CONDITION  FOR    RIGHT  TURN   IS   SHOWN 
ALL   RUNS  START   AT  GEOGRAPHIC   ORIGIN,  i.e., 

SHIP'S  COORDINATE  ORIGIN   IS  LOCATED   AT 
GEOGRAPHIC   COORDINATE  ORIGIN  AT   r*0.0. 
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The   equations   of  motion   can   then  be    approximated  by 


SURGE 

m(u  -  vr)  =  X 

SWAY 

m(v  +  ur)  =  Y 

YAW 

Iz  r      =  N 

Navigation 

(1) 

(2) 
(3) 


x     =  u   cos    ¥   -   v  sin   ¥  (4) 

y      =  u  sin   ¥   +  v  cos    ¥  (5) 

yo  

v     ="\(u2   +  v2  (6) 

3      =  tan"1    (2L-)  (7) 

In   the    above   equations: 

m  =  mass   of  the    rigid  ship 

I  =  moment   of   inertia  about   the    z-axis 
z 

u  =  v      cos    3    =  velocity   in   the  x-direction    (SURGE) 

s  ■* 

v  =  V     sin    3    =  velocity   in  the   y-direction     (SWAY) 
s 

r  —  $ -   angular  velocity   about   the    z-axis 

¥   =  heading    (yaw)    angle    (see   Fig.    1) 

3    =  drift   angle    (see   Fig.    1) 

5  =   rudder  angle,    +  5  produces    a  left   turn    (see   Fig.    1) 

Also: 

=   /xQdt  u  =j\idt  r  =jr 

=  yyodt  V=/vdt  f-/f 


xo 


Yo 


rdt 
fdt 
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The    right-hand-sides    (RES)    of   the   equations    of  motion 
have   been  designated  X,    Y,    N.      These    symbols    represent 
summations   of    forces    (X,Y)    and  moments    (N) ,    and  each   summation 
may   consist  of  many    terms.      For  the   initial   simplified  model 
we   have    discarded   all  terms   in   the   summations   except  the 
few  that    are   obviously  major   contributions.      The   terms    re- 
tained  are    considered   sufficient   to   produce    a  working 
simulation  which   should  provide   qualitative    agreement  with 
reference   data.      At   a   later  point   the   model  may  be    refined 
by   adding   terms    and   adjusting  terms    and   adjusting   coefficients 

After   reduction   of   the    right-hand-sides,    these    force 
summations   become: 

X  =    GT    cos    6    -    C      u  (8) 

Y   =    GT    sin    6    -    C      v|v|  (9) 

N   =   -(GT    sin    6)1      +   C_      vlvll  ,  ,n> 

o  Dy       '     '    u>  (10) 

where:       (see    Fig.    2    also) 

GT   is   the    gross   thrust   in  pounds 

5   =    5      =6      =   angle   of   all   thrust  vectors 
s  P 

C        =    coefficient   of   drag  in   the    direction   of   surge 

C_      =   coefficient   of   drag   in   the    direction  sway 

1       (see    Fig.    2)  =   thrust    lever   arm 
o  3 

1       (see    Fig.    2)=   sway-drag   lever   arm 

In    further  explanation   of  equations    8,    9,    10,    note 
that  we   have    chosen   to   model   the    drag   forces    as    lumped 
parameters    acting   at   two   points   P    and  P "    (see   Fig.    2)  . 
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P  -  P    ARE  EQUIVALENT   FORCE   CENTROIDS 
Figure  2.    FORCES  ACTING  FOR  A   RIGHT   TURN 
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The  length  1   is  the  distance  from  the  stern  to  the  center 

of  gravity,  and  1   is  the  component  distance  (in  the  x- 

direction)  from  the  center  of  gravity  to  the  points  P-P". 

1   is  a  fixed  parameter  for  our  initial  studies  (it  could  be 
o 

changed  by  ship  loading)  but  1   is  an  adjustable  distance 
which  we  will  use  as  a  parameter  to  make  the  model  fit  test 
data. 

The  coefficient  C_„  is  a  lumped  parameter  equivalent 

for  the  skin  friction  effect  (slender  body  theory) .   The  term 

2 

Cn„u  models  the  hull  drag,  aerodynamic  effects  and  hydro- 
dynamic  effects.   The  bubble  drag  is  not  included  in  the 
initial  model  (will  be  addded  later)  because  our  initial 
efforts  are  aimed  at  evaluation  of  our  hybrid  structure, 
and  a  simple  model  above   critical  speed  will  suffice  for 
this.   In  like  manner,  in  eqn.  9,  the  coefficient  Cn  models 
the  cross  flow  drag  effects.   Note  that  the  added  mass  effect 
has  been  temporarily  omitted  since  its  effect  is  also  of 
second  order. 

Substituting  equations  8,  9,  10  into  equations  1,  2,  3 
provides  the  basic  equations  of  motion  for  the  3  degree  of 


(GT  cos  5  -  C^u2)  +  vr  (11) 

(GT  sin  6  -  CD  v|v|)  -  ur  (12) 

-(GT  sin  6)1   +  (C_v|v|) 1   ; 
o      Dy  '  '   a)  ' 


Equations  11,  12,  13  are  applicable  to  any  SES  craft. 
Selection  of  parameter  values  to  be  used  in  the  equations 


freedom 

model 

• 

u 

= 

1 
m 

• 

V 

= 

1 
in 

• 

r 

= 

1 

I 
z 
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would  depend  on  the  particular  craft,  but  the  procedures 
used  to  identify  the  lumped  parameter  values  are  independent 
of  the  craft  type . 
3.   IDENTIFICATION  OF  CRAFT  PARAMETERS 

In  order  to  evaluate  the  parameters  needed  for  the  equa- 
tions, the  craft  must  be  available  for  measurement  and  test 
(or  suitable  data  must  be  available) .   The  parameters  m,  I  , 
1   are  constants  which  must  be  known,  and  GT,  S,  u,  v,  r 

can  be  measured  in  steady  state  (if  needed)  .   Then  C_w  1  , 

DX   go 

C_   can  be  calculated  as  is  shown  below: 
A.   Surge  Drag  Coefficient,  C 

The  craft  is  run  on  a  straight  course  with  GT  = 
constant  and  given  rudder  conditon  5 .   When  steady  state  is 
reached  u  is  measured.   Then 

fi  -  °  ■  -S-  (GT  cos  s  '  Sx^  (14) 

from  which 

CDX  =  GT  cos  5/u2ss  (15) 


where 


GT  =  gross  thrust  in  pounds 
u   =  steady  state  speed  in  feet/sec 
5  =  rudder  angle 
B.   Sway  Drag  Moment  Arm,  1 

Let  GT  =  constant  and  5  -  <5,  ,    a  fixed  angle.   When 
the  craft  reaches  steady  state  in  the  turn  the  sway  equation 
becomes 

v  =  0  =  -~(GT  sin  5  -  CD  v|v|)-  ur  (16) 


17 


from  which 

2 
C^  v  =  GT  sin  5  -  m  u  r  (17) 

Dy 

The  steady  state  yaw  moment  equation  becomes 

r  =  0  =  ~y~     -(GT  sin  5)IQ  +  (C  v|  v| )  1  1       (18) 

from  which 

C   v2  =  (GT  sin  5) (1o/1  )  (19) 

Combining  equations  17  and  19 : 

GT  sin  6(1-  1o/1  ) 

r  =  2-  (20) 

mu 

or  if  we  wish 

1   =  i   GT  sin  5 (21) 

oj     o   GT  sin  5  -  m  u  r  v  ' 

From  Fig.  2  and  eqn .  20  it  is  clear  that  the  turning 

rate,  r,  is  a  function  of  the  distance  1  ,  with  all  other 

parameter  values  fixed.   Since  it  is  customary  to  command 

a  turning  rate,  if  we  choose  r  we  can  then  evaluate  1  .  For 

the  XR-3,  using  steady  state  data  obtained  with  our  6  DOF' 

simulation, 

r  =  -0.0266(1  -  12/i  )  radians/sec  (22) 

u 
A  plot  of  this  relationsip  is  shown  on  Fig.  3. 

C.   Sway  Drag  Coefficient,  C 


Again  assuming  steady  state  conditions  with  fixed 
thrust  and  fixed  rudder,  the  sway  equation,  eqn.  17,  provides 

r        =  ~m  u  r  +  GT  s^n  -  f)-x\ 

LDy  '         2  U6] 

2  v 


r    (rod /sec) 
A 


'.]%)   .3   - 


4°4)    .2  - 


7%)  .1    - 


Figure  3.  PLOT  OF  RELATIONSHIP 
BETWEEN  r   AND      ! 


r  =  -.00266    (1-    /,    ) 
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Figure  4-A.    FLAT  TURN  DYNAMICS   OF   THE   XR-3,  PLOT 

OF   VELOCITY  vs.  TIME. 

r=    .33   RAD/SEC    ( L  =  I  ft) 
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Figure  5-A.    FLAT  TURN  DYNAMICS  OF  THE   XR-3  PLOT 
OF  VELOCITY  vs.  TIME. 
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Figure  5-B.   FLAT  TURN  DYNAMICS  OF  THE  XR-3   PLOT 
OF  X  vs.  Y. 

r=  0.1  RAD/SEC    (  I     =2.5  ft) 
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For  a  given  turning  rate,  r,  we  know  all  values  except  v, 
which  can  be  calculated  using  the  surge  equation,  eqn.  11, 

to  provide: 

C-^u   -  GT  cos  5 

v  =  — (24) 

m  r 

4.   DIGITAL  SIMULATION  OF  THE  SIMPLIFIED  XR-3  EQUATIONS  IN 

3  DEGREES  OF  FREEDOM 

The  equations  of  motions  as  previously  derived  were 
programmed  into  the  IBM  360/6  7  computer.   A  fixed  step  size 
was  used  to  verify  the  program  and  validate  the  model. 
Using  a  turn  rate,  r  =  .33  rad/sec,  and  an  integration  step 
of  0.1  seconds  the  simulation  was  run  for  120  seconds  of 
problem  time,  then  repeated  for  an  integration  step  of  0.01 
seconds.   Comparison  of  results  showed  that  an  integration 
step  of  0.1  seconds  was  suitable,  all  variables  in  the  two 
runs  agreeing  to  within  1% . 

Figure  4A  and  4B  show  plots  of  u  vs  time  and  x  vs  y 
for  a  flat  turn  with  r  =  .33  rad/sec  (1   =  1  ft.) .   Figures 
5A  and  5B  show  similar  plots  for  r  =  0.1  rad/sec  (1   =  2.5  ft 
A  number  of  additional  runs  were  made  but  results  are  not 
given  here.   Based  on  these  runs  we  feel  that  the  program 
has  been  satisfactorily  verified,  i.e.,  it  is  solving  the 
equations  correctly. 

The  following  observations  seem  appropriate: 

1.  Eqn.  22  and  Fig.  3  are  at  least  qualitatively 
correct. 

2.  The  effect  of  added  mass  is  neglegible. 
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3.  We  have  gained  confidence  that,  given  the  meas- 
ured perfomance  of  any  SES  in  an  approximate  flat  turn,  and 
given  values  for  the  parameters  GT ,  5,  u,  r,  m  and  1  ,  we 
can  determine  a  value  of  1  which  provides  a  simulation 

OS  c 

matching  the  measured  performance  reasonably  closely. 

4.  If  the  simulation  is  performed  with  an  initial 
step  in  the  thrust  angle,  5,  approximate  values  for  the  system 
time  constants  can  be  evaluated  from  the  resulting  data. 

For  the  XR-3  operating  in  a  flat  turn  the  time  constants  are: 

Variable  Time  Constant 

u  11.0  sec 

v  0.8  sec 

r  1.2  sec 

6  0.9  sec 

To  validate  the  3  DOF  model,  additional  runs  were  made 
using  the  digital  6  DOFOceanics  program  as  modified  by  our 
NPS  group  for  the  XR-3.   Also  the  3  DOF  model  was  run  with 
3K-SES  data,  and  the  5  DOF  model  of  the  3K-SES  was  run.   For 
the  XR-3  tests  a  step  of  rudder  of  0.5  radians  was  suddenly 
applied  with  craft  speed  at  20  knots.   For  the  3K-SES  tests 
and  11   rudder  was  applied  by  ramping  up  for  1.0  second  and 
holding  at  the  11   value.   Craft  speed  was  60  knots.   Com- 
arison  of  the  3  DOF,  6  DOF  and  3  DOF  Hybrid  runs  for  the 
XR-3  is  shown  on  Fig.  6.   Correlation  is  not  as  bad  as 
appears  by  inspection.   Turning  diameter  and  time  for  a  360 
turn  are  almost  indentical.   We  presently  attribute  the 
differences  to  the  fact  that  the  3  DOF  model  does  not  include 
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Figure  6.    TURNING   CHARACTERISTICS   OF  THE    XR-3    RESPONSE 
TO  LEFT   RUDDER,  0.5   RADIAN  STEP,    20   KNOTS. 
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a  term  for  bubble  drag. 

Figure  7  compares  the  3  DOF  run  with  the  5  DOF  run 
for  the  3K-SES.   The  results  are  almost  identical,  despite 
simplifications  in  the  3  DOF.   This  gives  confidence  that 
our  modelling  technique  will  work  well  when  we  hybridize. 

5.   HYBRID  SIMULATION  OF  THE  3  DEGREE  OF  FREEDOM  MODEL  OF 
THE  XR-3 


In  developing  a  hybrid  model  various  arrangements  and 
structures  (programming)  are  possible.   The  philosophy 
chosen  for  our  initial  hybrid  model  is  shown  in  the  block- 
flow  diagram  of  Fig.  8.   Details  of  the  programming  are 
given  in  Appendix  A  for  the  simplified  XR-3  model. 

The  simplified  XR-3  model  was  installed  in  the  hybrid 
system  and  tested.   Comparison  of  results  with  those  ob- 
tained by  digital  simulation  showed  almost  identical  per- 
formance, differences  being  less  than  2%. 

In  developing  the  hybrid  simulation  a  major  concern 
has  been  the  amount  of  time  delay  introduced  into  the  loop 
by  the  A/D  and  D/A  conversions  and  by  the  time  required  to 
perform  the  digital  computations.   We  have  monitored  these 
delays  in  the  3  DOF  model  with  results  as  follows: 

1.  The  fundamental  loop  delay  using  D/A  for  u,  v,  r, 
x,  y  and  A/D  for  u,  v,  r  was  33  ms . 

2.  A  lower  language  subroutine  exists,  called  ADDA, 
which  can  be  used  to  replace  the  Fortran  subroutines 
for  A/D  and  D/A.   This  can  reduce  the  delay  by  25%. 
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3.   Additional  delays  are  introduced  by  using  A/D 

for  thrust  (GT)  and  thrust  angle  (6)  f76%  increase] 
and  by  use  of  the  Graphics  Signal  (Gated)  conversion 
of  the  variables  [6  7%  increase].   Total  loop  trans- 
port delay  using  D/A,  both  A/D  units  and  the  Graphics 
signal  conversion  was  approximately  7  8  ms. 
We  know  how  to  reduce  the  loop  delay  substantially,  but 
in  view  of  the  system  response  time  this  may  not  be  necessary. 
One  possible  source  of  difficulty  in  the  hybrid  simulation 
is  the  cross  coupling  between  the  equations.   When  operating 
near  steady  state  conditions  during  small  yaw  angle  maneuvers 
the  linearized  equations  indicate  that  the  roots  of  the 
characteristic  equation  may  be  very  lightly  damped. 

Consider  straight- ahead  motion  (thrust)  at  constant 
speed.   The  linear  model  sway  and  yaw  equations  are: 

* 

2     22  2     o  3 
x3  =  a32x2 


where 


a  r        It 

x  =   Au,  Av,  Arl  ,  and  Au  = 


a22  "  fl(x2} 

a32  =  f2(x2) 
Combining  these  equations,  the  characteristic  equation 

of  the  simplified  3  DOF  model  is 

2 
s   +  anos  +  u  a__  =  0 
22      o  32 

The  locations  of  the  roots  of  this  equation  depend  on  the 
values  of  a^-  and  a.,- ,  both  of  which  are  functions  of  the 
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sway  velocity.   For  small  values  of  sway  velocity  the 
roots  are  complex  conjugates  and  approach  the  origin  of 
the  s-plane  as  the  sway  velocity  approaches  zero.   For 
larger  values  of  sway  velocity  such  as  those  encountered 
in  a  turn  maneuver  the  roots  are  real  and  negative,  pro- 
viding a  well  damped  response  for  yaw  and  sway  motions. 
This  appears  to  be  similar  to  the  Dutch  Roll  effect  in 
aircraft. 

For  straight  ahead  runs  using  the  digital  program  for 
the  3  DOF  model  and  setting  the  initial  value  of  v  to  zero, 
no  sway-yaw  motions  are  observed.   However,  if  we  set  a 
small  initial  condition  on  v  a  transient  oscillation  is 
observed  in  both  sway  and  yaw  as  shown  on  Fig.  9.   This 
verifies  the  existence  of  the  phenomenon,  and  we  expect 
to  encounter  it  in  the  hybrid  since  the  equations  are  the 
same.   However,  in  the  hybrid  there  is  a  possibility  that 
the  roots  may  move  into  the  right  half  s-plane  due  to  either: 

1.  DC  unbalance  in  the  analog  computer  integrators,  or, 

2.  Excessive  loop  delay  from  A/D,  D/A  and  digital 
computations . 

The  first  of  these  possible  causes  can  be  avoided  by 
careful  balancing  of  the  operational  amplifiers  before 
running.   The  effect  of  the  loop  delays  can  be  controlled 
somewhat  by  minimizing  the  number  of  terms  requiring  digital 
computation. 
6.   STATUS  OF  THE  3  DOF  MODEL  STUDIES 

These  appear  to  be  both  complete  and  successful.   Both 
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qualitative  and  quantitative  agreement  have  been  achieved 
with  the  5  DOF  digital  simulation  as  a  reference.   Future 
studies  will  be  aimed  primarily  at  the  5  DOF  model  of  the 
3K-SES  except  when  we  wish  to  correlate  computer  data  with 
specially  designed  experiments. 
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IIB.   THE  5  DEGREE  OF  FREEDOM  MODEL 

1.  INTRODUCTION 

Having  developed  a  simplified  3  DOF  model  of  the  XR-3, 
we  now  proceed  to  use  the  knowledge  thus  obtained  and  apply 
it  to  the  3K-SES.   We  do  not  plan  to  develop  a  5  DOF 
hybrid  model  of  the  XR-3.   Our  procedure  is  to  install 
the  NSRDC  5  DOF  program  for  the  3K-SES  in  our  IBM  360/6  7 
computer,  generate  reference  curves  with  it,  obtain  any 
numbers  we  need  for  our  5  DOF  hybrid  simulation,  then 
verify  and  validate  the  hybrid  model  of  the  3K-SES. 

2.  THE    DIGITAL    SIMULATION   OF   THE    3K-SES    IN    5    DEGREES    OF 
FREEDOM 

This  program  was  obtained  from  NSRDC.   It  has  been 
written  for  a  CDC  machine  and  had  to  be  converted  to  the 
IBM  version  of  Fortran  for  our  computer.   This  has  been 
completed,  the  program  is  operational  and  as  far  as  we 
know  has  no  obvious  "bugs".  We    are  familiarizing  ourselves 
with  the  program,  and  exercising  it  to  obtain  data  we  will 
need  in  developing  the  5  DOF  hybrid  model,  and  for  validating 
it. 

Details  of  the  conversion-problems  encountered,  their 
solutions,  etc.  are  discussed  in  Appendix  B. 

3.  HYBRID    SIMULATION    OF    THE    3K-SES 

This  is  under  development  at  the  present  time.   Hope- 
fully we  need  only  insert  a  new  set  of  constants  in  the 

34 


equations   of  the   hybrid  XR-3  model,    but   undoubtedly  other 
adjustments    and  modifications  will  be   needed,   parti cularily 
in  the   trade-off  of  digital   vs    analog   computation   in   order 
to  perform  a   real-time    simulation. 

4 .      STATUS 

Progress  has  been  satisfactory.   It  is  too  early  to 
predict  future  developments. 
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III.       SUMMARY    AND    CONCLUSIONS 
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The  NPS  Hybrid  Computer  consisting  of  the  XDS-9  300 
digital  and  the  CI  5000  analog  computers  has  been  pro- 
grammed to  simulate  the  3  DOF  Flat  Turn  maneuvering 
dynamics  of  the  CAB  type  Surface  Effect  Ship. 

The  XR-3  has  been  programmed  with  simplified  modelling 
of  the  drag  forces  and  thrust-input  force  on  the  hybrid 
computer  for  simple  flat  turn  maneuvers.   An  evaluation 
of  the  hybrid  simulation  was  carried  out  by  comparison  of 
the  time  histories  of  the  velocities  and  the  navigation 
variables  against  those  values  obtained  from  the  3  DOF 
digital  simulation  for  the  same  maneuvers.   In  addition, 
the  XR-3  hybrid  3  DOF  results  were  compared  to  the  6  DOF 
loads  and  motion  output . 

The  result  of  this  comparison  was  that  the  hybrid 
computer  simulation  produced  a  time  response  to  a  step 
rudder  input  that  produced  a  navigational  turning  plot  that 
was  in-between  the  all  digital  3  DOF  and  all  digital  6  DOF 
trajectories. 

Additionally  the  3  DOF  simplified  RHS  (right-hand- side) 
digital  program  was  operated  for  the  same  turning  conditions 
that  was  used  for  the  5  DOF,  DB.SIM5D   Program  and  the 
comparison  indicates  close  agreement  in  the  navigational 
trajectories  for  an  11  degree  vectored  thrust  angle  full- 
turn  maneuver. 
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This    feasibility   study  has    shown   that   the   hybrid 
computer  implementation    for   the    5    DOF   3K-SES    in   real- 
time   is    possible   with   the    available   equipment    at   the   NPS 
Computer  Laboratory.      Real-time    computation   is   of   course 
only   approximated  when   using  the   hybrid   computer.      The 
loop   transport   delay   inherent   in  this    simulation   is    due 
to   the    signal   conversion    (A  to   D   and  D   to   A)    and  the 
necessary   digital    computation    for  the    acceleration   forces 
and  moments.      The   magnitude   of  this    delay  time   has    a   fixed 
minimum  value    as   a   result   of   signal   conversion   and   a  maximum 
value    due   to   righ-hand-side    calculation   required   for   the 
acceleration   added  to   the    fixed  time    due   to   the   signal    con- 
version.     Real-time    computation   accuracy   is    determined 
principally  by   the  magnitude   of   this   maximum  value    for   the 
loop   delay. 

The    3K-SES   was   approximated  with   the    simplified    3    DOF 
Flat   Turn   equations    and  the    results  were    compared.       It 
was   shown   that   the    full  maneuvering   trajectories    for   an 
11   degree    thrust   angle   compared  very   closely  with   the   5    DOF, 
DBSIM5D  Program  results.      These    results    indicate    that  only 
a   few  terms   on   the    RHS   may  be    all    that   are    required   to    simu- 
late   the    3K-SES    for  the   type   of  maneuvering   control    studies 
to  be    conducted   in  the   future   on   the   hybrid   computer. 

The    conclusion   of  this    feasibility   study   is    that    suc- 
cessful implementation  of   a  hybrid   computer,    real-time 
simulation   of   the    DBSIM5D   Program  has    a  very  high   prob- 
ability.     The    degree   of    complexity   necessary    for   the    RHS 
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computation  will  be  the  limiting  factor  on  the  loop  delay 
that  controls  the  real-time  accuracy. 
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APPENDIX  A.   HYBRID  COMPUTER  PROGRAM  FOR  THE 
3  DOF  (FLAT)  TURN  EQUATIONS 

A  general  schematic  diagram  of  the  hybrid  structure 
used  is  shown  in  Fig.  A-l.   The  listing  of  the  digital 
program  for  the  XDS-9  300  computer  is  also  included  in  this 
appendix,  and  the  analog  computer  connections  are  shown 
on  Fig.  A- 2. 

As  may  be  seen  from  the  listing  of  the  digital  pro- 
gram, the  right-hand-sides  (RHS)  of  the  acceleration 
equations  for  the  3  DOF  are  computed  on  the  XDS-9300, 
which  also  does  the  magnitude  scaling  for  the  analog 
computer  (CT-50  0)  voltages.   In  addition  the  XDS-9  300 
provides  the  required  control  logic,  timing  and  temporary 
storage  of  computed  variables. 

The  graphics  terminal  (AGT-10)  provides  real  time 
display  of  the  desired  output  and  control  variables,  as  well 
as  their  maximum  and  minimum  values.   Piloting  of  the  craft 
is  carried  out  at  the  graphics  terminal.   Additional 
permanent  record  outputs  of  navigation  x  and  y  are  obtained 
with  an  XY-pen  recorder,  and  plots  of  other  system  variables 
are  obtained  on  an  8  channel  Brush  recorder. 

A  sample  of  the  x-y  output  for  the  3  DOF  simplified 
XR-3  craft  is  shown  on  Fig.  A-3.   The  plot  starts  with  a 
right  turn  (5  =  -0.5  radians)  at  t  =  0 ,  using  a  terminated 
ramp  input.   The  rudder  was  held  at  5  =  -0.5  for  two  turns, 
then  reversed  to  5  =  +0.5  for  a  half  turn,  followed  by 
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Figure  A-l.     BLOCK  DIAGRAM  OF  THE  HYBRID  STRUCTURE. 
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6  =  0.0,  then  again  <5  =  +0.5  for  a  quarter  turn,  and  finally 
a  right  turn  using  6  =  -.29  radian  for  two  complete  turns. 
The  total  time  for  these  maneuvers  with  the  XR-3  was  approxi- 
mately  3-j  minutes. 

Note  that  the  CI-500Q  is  used  only  for  analog  inte- 
gration of  the  accelerations  to  obtain  the  desired  velocities 
for  the  digital  computation.   Displacement  is  computed  in 
the  digital  computer  using  the  trapezoidal  rule  for  integratior 
as  done  on  the  Oceanic' s  6  DOF  Loads  and  Motions  Program. 

Note  also  that  the  present  structure  and  programming 
use  only  existing  hardware  and  we  consider  them  a  feasi- 
bility study.   If  results  show  that  we  have  been  successful, 
then  we  must  prepare  for  manual  control  studies,  and  this 
may  require  development  of  a  more  realistic  "cockpit"  to 
simulate  conditions  in  the  control  room  of  the  3K-SES . 
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Figure  A-3.   SAMPLE   OF  THE  HYBRID 
OUTPUT  x-y   PLOT  FOR 

THE  3  D.O.F.  XR-3 
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APPENDIX  B.   INSTALLATION  OF  THE  3K-SES,  5  DOF 
PROGRAM  IN  THE  IBM  360/6  7  COMPUTER 

In  the  fall  of  19  7  8  the  3K-SES,  5  DOF  data  based 
program  called  DBSIM5D  was  delivered  to  the  Naval  Post- 
graduate School,  at  which  time  a  commitment  was  made  by 
the  NPS  group  to  attempt  a  hybrid  computer  simulation  based 
on  this  program. 

In  the  spring  of  19  79,  DTNSRDC  shipped  the  card  deck 
and  the  output  data  for  a  typical  3K-SES  run,  obtained  by 
DTNSRDC  on  a  CDC  computer.   The  deck  was  to  be  used  in 
installing  the  program  in  the  IBM  360/6  7,  and  the  output 
data  provided  information  for  verification  and  validation. 
The  major  problem  in  the  process  of  installation  was  the 
conversion  from  CDC  Fortran  to  IBM  Fortran.   With  the  help 
of  Mr.  Roger  Hilleary  of  the  NPS  Computer  Facility,  the 
interactive  graphics  terminal  was  used  to  edit  the  program. 
After  some  50  man  hours  of  effort  the  first  successful 
batch  run  was  made  in  April  19  79. 

Complete  success  in  the  conversion  did  not  occur  until 
a  correction  was  added  to  the  newly  programmed  TIMER,  to 
allow  for  a  stop  at  the  desired  stop  time  rather  than  a 
premature  stop  time  due  to  the  timer  program. 

To  validate  the  IBM  version  of  the  program  we  dupli- 
cated the  sample  run  for  the  3K-SES.   This  run,  as  provided 
by  DTNSRDC,  was  for  a  left  turn,  using  an  11  degree  effector 
angle  and  a  speed  of  60  knots.   The  sample  run  terminated 
after  14.5  seconds  of  real  (problem)  time.   The  11  degree 
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effector  angle  was  programmed  as  a  ramp  which  terminated 
after  1  second,  and  remained  constant  at  11  degrees  for 
an  additional  50  seconds. 

A  comparison  of  the  output  of  the  IBM  360/6  7  version 
with  the  sample  run  data  from  the  CDC  computer  showed  that 
the  recorded  variables  were  in  complete  agreement  within 
the  accuracy  differences  of  the  two  computer  word  length. 

The  problems  of  conversion  were  simple  but  time 
consuming  and  may  be  listed  as  follows: 

1.  6  character  variable  name  limit  on  the  IBM  com- 
puter as  compared  with  a  7  character  limit  on  the 
CDC. 

2.  Single  statement  limit  on  the  IBM  computer  versus 
multiple  statement  on  a  single  card  for  the  CDC. 

3.  Identification  of  temporary  files  used  for  the 
output  data. 

4.  Identification  of  the  timing  control  and  pro- 
gramming the  IBM  computer  timer  to  provide  the 
functions  required  for  operation  of  the  program. 

5.  Addition  of  the  proper  job  control  language  (JCL) 
cards  to  compile,  link  and  compute. 
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DAT 


AreRTRA\ 
1 

3 

4 

5 
6 

7 

3 

Q 

to 

11 
1? 
13 
14 

15 
16 
17 

IS 

19 
20 
21 
22 
23 
24 
25 
26 
27 
2  H 
29 
30 
31 
32 
33 
34 
35 
36 
37 
33 
3? 
40 
M 
42 
43 
<+* 
45 
**6 
47 
4? 
if? 
50 
51 


CI 
C/I 
CIL 
CIV 

EQ 
C(A 
C(A 
C<A 
C(  4 

CA 

ca 

CA 

C* 
CA 
CA 
CA 
CA 

EN 
p» 

CA 
EN 
F5 
CA 
1  = 
CA 
IV 
IV 

IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
CA 

IP 

IC 

i: 

\D 

I  N' 

x  = 

TI 
A  A 
P  = 

0  = 
PH 


VrZNS 

L(24 
A4(2 
IEwB 
U  I  V  A 
A(2) 
A(7) 
A  (  1  2 
D(7) 
LL  0 
LI  S 
PO  15 


LL  R 
LL  C 
LL  0 
LL  0 
LL  D 
LL  0 
C^CE 
R^AT 
LL  T 
C?DE 
R^AT 
LL  T 
i  0  0  0 
lL  0 
IEwB 
IEwB 
IEwB 
IEWB 

IEWB 

T  r ,.,  a 

IEaB 
IEwB 
IEWB 

IEaB 


LL 
R=C 

9A  =  0 
5  =  0 
A»12 

0  =  10 

=  0 

-ICO. 

*E*0.< 

AA  =  0  . 

CO 

CO 

e  *  2 1  o 


I9N  I 

)  ,TIT 

4),IL 

(20) 

LENCE 

#321) 

#861) 

>  /  B 1 1 

#A71) 

AU<  .0 

,  .0) 

ESET( 

3  m,  o  \j  j 

G  T  N  I T 
G  !  v  I  T 
GIN  IT 
TINIT 
(96/1 

Exre( 

(96*2 
(  » 
EXTS( 

000 
ELAY 

( 1 )  =  1 

(2)»I 
(3)«I 

(  4  )  =  I 
(5)  =1 
(6)»I 
(7)«I 

(3)»I 
(9)=I 
(10)  = 
(  11  )■ 
(12)  = 


TrxT(29) , ILAP(2 
LE(24),AC( 1G)/A 
A5(24), I^A6(24) 

(AD(  1)/  ah  )  ,  (A 
i  ( AA(3),  331 ) >  ( A 
i ( AA( 3) #571 ) t   ( A 

1  )  ,  (  A  Q  (  4  )  ,  A  4  1  )  / 

/  (  AO(  8)/  A81  )  /  (  A 
/•0j.0/.0/»0/»0 
(  4HP001/  .C^-D3 

1000) 

(2i IVIEWA/2C* IE 
(2/  IVIEWE/20/ IE 

(2# i°LeT/ ioj :ep 

( li  ITEXT,29,  TER 
/TITLEC) 

REAL    TIV 
l/TrTLE0#24,6il 

#  TITLED 
MeTlBN    ANALYSI 

1/TITLE1,24,9,1 


4)iP0T3(24 ),°5TK(24), IM 
A(12)  , I  LA  1 (24) , ILA2(24) 
,  lLA7(24), ILAB(24)/ IPL9 

D<2)/ A21 ) , (A0< 3)/A31 )i ( 
AC*  ),3'+l  )j  {  AA(5)  ,  B5l  )  ,  < 
A (9)/331> j ( aa( 10)/391 )/ 
(AD<5)/ A51 )* < AD(6)/A61) 
D(9)/A91  )/  (  A0(  lO/AiGl) 
/  .0/  .0* ,0j «0/ t0# .0) 
03/  .0/VHPC05/  «0#  4HP007# 


P) 
R) 
) 
) 


A3E(30 
>  I„A3( 
T(IO)*- 

AA(  1  )/ 
A  A  (  6  )  » 
(  A  A  (  1 1 


•J  sj    I 

24), 
iVZEMACa 

311)* 

355  )/ 

) ,  3 1 0 1 )  J 


•Qj 4HP011/ «0i 


E 

/  J/  -3  /    1 


^•^•TERj 


S    PRQGRA^ 
,3/3, IER) 


^E*D< 
PAC'<( 
PAC<( 
PACK( 
PAC'<( 
PAC'<( 
PAC'<( 
PAC<( 
PAC<( 
Tdac'< 
IPAC'< 


1/10) 

0.0/1.0*0) 
0.0/0.5/ 1) 

.5/C  0/1 


0.5/0.0/0; 
1.0*0*0/ 1) 
0.75/1 .0/0 

0.73/ 0.5/1 
(0.5/3.75/ 
(1.0,0.75/ 


) 
) 
0) 

1  ) 


2/  IV IEWB/ 12/1 j IER) 
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date 


G=32.2 

M  s  1  5  C  • 
C^X=.3 
C0Y=1555.0 
LW*1  .0 

J=l 

IZ=5333. 
A3  2  =  L0 
All*  #  <+5 

AUD3T««45 

AV=.Q2C 

AVD?ts.C2C 

AP-O.C100 

aPD°T=0.  010C 
AXC=25. 
AY0  =  25  . 
AT =3. 75 


TX=0«C 
TY=-6000t 
fc5c    C*LL    DTINIT<1*ITEXT/29*IER) 

501    F?RN'AT(  'PRESENT'  i2F10.2*  ' 

C  '  •  /  F  6 . 3 ) 

525    F^R^ATf  i.yAXl^w'!*'  #2FXC  »2j  ' 

C  SF6.3) 

5  26    F  S  R  w  A  T  (  •  v  J  m  t  v  >j  m  i  ,  ^  r  j_  g  ,  p  ,  i 

C  *  iF6.3) 

5C7    F9RMAT(  'ORESE^T»/F8.2/FlQ.3j 

CF5  .2i  •  '  ,F5.3) 

527    F?P^at(  'MAXImUm'jF8«2<»F1Q»3j 

CF5.3; •  > »F5.3) 

5  23    F^PviT;  »  MiNIf-UMl*r  8«2*FlO»3i 

CF5 »3i '         '  *F5»3) 
511  F9RMAT< 'PRESENT' •F7»3jF10«3j 

C     ' #F6.3) 

529  FSRMAT< 'MAXIMUM' /F7» 3* F10»3j 
C     *iP-6.3) 

530  F9RyaT(  »!-!IM^Um,#F7.3*F10»3j 
C     SFA.3) 


516  F?P^AT<F6 ,£• 


';-^.?> ' 


7»2i ' /' jF5»1 


F5.2. 


F^.p,  i 


*F4-2j ' 


i  .  r 


•  c 


/F3«2^  * 


f  -  6  •  cf 


/ '  '/F8.3 


i 


1,-4.2,'  >  ,  -r3.3;  ' 


*  F 1 C .  2  #  ' 


• #F1C«2* ' 


'   ' >r5.Pi ' 


•  j  F  8  •  3  ,  • 

'  #F8#3j ' 
'  jF8.»3#  ' 


'/F5i2) 

pasirieM' ) 


ENCeCE(96j5CC/ I _A1 ) 

SCO  F?R"AT{ ' 

CALL  TEXTCdi  lLAl#2*#3#lil#3j  IER) 

IP(  IcR.NF.O)?w:r  =  'jT(  1C1  )  IERj  '     P6SITI9N' 

E\C=CE(56i502/ ILA2) 
502  F?P^at; '  X  Y  ZZ         PSI 

CDH I  Tur  TA '  ) 

C^LL  TEXie(  l/lLA2/2^^>l>l/3/IE3) 

ENCeDE(96/505f ILA3) 
505  F9RMAT( '  VELOCITY') 


F6.3, 


,F6«3j ' 


F6«3i 


1  *  ra  •  3. 
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DATE 


509 


512 


11 


CALL  T 

ENCODE 
17  F9R^AT 

CP 
CALL  T 
I F  (  I E  ^ 
ENC50E 
F?R^AT 
CALL  T 
IF( IER 
ENC8DE 
F9RMAT 

C9T 
CALL  T 
in  IER 
ENCODE 

14  FORMAT 
C  I G  A  T  I  Q 

CALL  T 
IF(  IER 
ENC6DE 

15  R9RN"A'r 

CaTE 
CALL  T 
IF( IER 
CALL  ••• 
CALL  S 
C8NTIN 
CALL  A 
CALL  R 
CALL  W 
X1=PSI 
I X  -  v  i 

AIX=X1 

I  f  (  a  i  y 

HEACsA 

A  I  X  1  s  1 
hEAO=A 
C9NTIN 
ARATE= 
AR9LL* 
A^ITCU 
S3EED= 
I F  ( H 1 . 

M 1  »  X  0 

IFCH2. 

H2sYC 

IF(H3- 

h3  =  ZZ 

IF(M4. 

H^=°SI 


R 


EXTg(1,!lA3/2^/11/1j1/3,IZR) 
.NE-O)SUT^UT(lOl)     IER/     '         VELOCITY' 
(96/517/ ILA^O 
(  '  u  V  W 

Q'  ) 
EXT9( 1/ I  LA 4* 2*/ 12* 1/1 j 3/ IER) 

.NE.0)8UTPUT(101)  IER*  '     .1ST  vAR  VEL3CITY' 
(96/509/ ILA5) 

(  '  ACCELERATION'  ) 

EXT6( 1/ILA5/24/19/1/1/3/IER) 
.MEtO)9UT°UT(101)  IER/  '    ACCEL  TITLE' 
(96/512/ ILA6) 
(  '  JD0T        VCrT         WD3T 

3D6T' ) 
EXTS( 1/ ILA6/24/20/1/ 1/3/ IER) 


RD  =  T 


VA,: 


TI*E 


CI 


700 
7C1 
702 


.NE.0)9UT°JT< 101)  IER/  • 

(96/514/ ILA7) 

(  •  C9NT?9L  INPUTS 

N  OATA' ) 

EXTS(l/lLA7/2^/27/l/l/3/ IER) 

.NE»0)e'JTPUT(101)  IER/  '    C3NTR9L  HEADING' 

(96/515. ILA3) 

('THRUST    RUDDER       OELT    ELAPSED 

DRIFT    ??LL    3=EED' ) 

EXT3(  lj  lLA8/2*/2S/  1/  1.3/  IER  ) 

.NE«0)8UTpUT(101)  IER/  •     3U3  u*3EL  C9NTR9L' 

RITECL9CK(Q) 

TARTCL^C'< 

UE 

DDA(AD/NAD/AA,NDA) 

EADCL5C*<(\) 

RlTECLdC'<(C) 

/6.  28 3 185 

-IX 

.LT.COG9  T3  ic 

IX*360. 

11 

•  +  AIX 

1X1*360. 

UE 

( 360»*R ) /6. 233135 

(360»*PHI )/6. 283185 

=(360»»THETA) /6. 233135 

L*. 597213 


3T.YG)G9  T9  701 


=  ITCH 


HEADING, 


3T.ZZ> 


'9  702 


GT.PSDG9    T9    703 
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DAT: 


7C3 
704 

7C5 
7  06 
707 
708 
7  09 
710 
711 
712 
713 
714 
715 
716 
71* 
713 
72C 
721 
722 
723 
724 
7?5 
726 
727 
728 
729 
73C 


IF(H5«GT#PHI )G3  TS  704 

H5*FHI 

IF(H6»GTtTHETA)G9  TO  70  5 

H6*T'HETA 

IF(H7.GT»U)G5  T5  706 

H7s'J 

IF( J8.3T.V)0e  T9  707 

N8*V 

IF{H9»GT»W)G9  T8  70S 

H9«» 

IF(H10«GT 

H 1 Q  *  R 

IF(W11.3T 


r-ll-° 

IF(hl2.GT 

hi2=: 

JF(H13.3T 
H13-UD9T 
IF(H14.GT 
W14*VCST 

IF(^15.3T 
n!5  =  ',DeT 
IF(H16«GT 
Hl6=RDeT 

IF(M17#GT 

Ml 7=°0^T 

IF(rll8-GT 
H18-GD3T 
I  =•  ( i- 3 C  -  LT 

~3C=xC 

Ic(^3! »UT 

H31=YC 

I'('-32.LT 

M32=ZZ 

IF  C-I33.LT 

H33aPSI 

IF(--3^.!_T 

IF(^35.LT 

H35=TH£TA 

[F(H36.LT 

H36*U 

IFM37.LT 

h37  =  V 

IF(M38.LT 

h38«* 

IF(-3?.LT 

H39sR 

IF(w40»UT 

H4QsP 

I  =•  ( :-  -+ 1  •  L  T 

H41=^i 

IF(H*2.LT 


R)G°  Tft  70S 

p ) g  e  t  e  7i c 

Q)G9  T8  711 
JDST)3?  TS  712 
VD9T)G9  Te  713 
W_09T)G9  T3  714 
RD9T)G8  T9  715 
POST). 39  TS  716 
00c"!33  T3  713 
XO  )  3^  TR  71° 
Y0)39  T9  720 
11)??    T3  721 
PS!)3?  T9  722 
PHI )Q3  T9  723 
THETAJG9  T?  72^ 
U)Ge  T?  725 
V)GO  "0  726 
w)G9  T9  727 
R)39  T?  72 S 


P)3?  T?  729 


)G9  T9  730 
D9T)3S  T5  731 


51 


DATE 


731 
722 
722 
73<+ 
725 
736 

651 


7  QO 

7?: 

732 


792 

795 

993 


H^2*uC9T 

lc  (H^3.L,T.VD9T)G8  T3  732 

H43=VD9T 

IF(M44.|„T»WDST)Ga  T3  "?33 

I^(H45.LT.RD?T)3e  T5  734 

H45=R03T 

IF(H46.LT»PD9T)38  T5  735 

p46=DD?T 

IF (Hf7.UT»3D9T)G9  Tr  736 

H47sQD8T 

C9NTINUE 

DELT*(DELT-MNi/60t  )  >/2. 

TI  -1F  =  7T^E  +  DELT 

Zs-C#5+A9I 

Ir( Z.GT. .5) 3?  T9  790 


IF(Z«LT»«»«5>G9 


731 


IF  (Z.GT.  .005)23  T3  792 
IF(Z.LT.-.005)G9  T9  792 

z=co 

29  T3 


Z=»5 


792 
792 


G3  Tg  732 

T=A'»A31 *1C0. 

1^(1.31.375. )G3  T3  793 

IF<T.LT.0«0)G9  T3  79* 

G  5  T  3  795 

T=275. 

39  T8  7Q5 

t  ■  o  •  a 

09  T9  795 

CONTINUE 

CONTINUE 

PAAP=0Q.2G*SIM3SI  ) 

PHI*R 

T'-tE'4=(jO?T 

*sPAAP 

P  =  PA  \P 

la  B  r  A  A  " 

WD9T=PAAP 

P09'=PAAP 

QD8T*PAAP 

U»AH*X*AU 

V*A2i*X*AV 

R=A3l*X*AR 

LD9T=((-1.5*(CDX/18C.)*j*J)  +  ((T/1S3.)*C^3(Z)5W*R 

V03T=( (T/M)*SIN(Z) >  -  ( (CDY/M)*V*A8S( V) )-J*R 

RD0T»(  (CDY/IZ)*LW*V*-ABS(V)  )-(  (T/IZ)»L9*3IN(Z)  )  +  (A22*j*V«(lPW/IZ)  ) 

BETAsATAN(V/'J  ) 

V3=  (  (U*U)+(V*V)  )**»5 
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DATE 


1  r*p 


1  ri(= 


653 
600 

107 

108 

109 
11C 

111 


6C1 
112 
113 


PSI=PSI+DELT*P 
X0D0T=U*C0S (PS 
Y0D5T=U*S  IN  (PS 
X0«X0+0ELT»XC0 
YO=YO+DELT*Y00 
311=UC9T/(AUD9 
B21»VC9T/ (  AV0A 
B31»RC9T/<ARDe 
371=x0/( ax0*1C 
BSlsY0/( AYO*10 
IF(SEkSE  SWITC 
CALL  W5L0 
CALL  ST9PCL9CK 
wPITE<  6i  100)  LP 
F9RMAT(9F9»4) 
CALL  STARTCL9C 
CALL  CgMPUTE 
C  9  N  7 1  n  'j  E 

ip<  noi.GT.-o. 
rce=iC3+i 
ir( ico.eo 

IP(  ICQ.EO 
Ir( ICO. EG 

IF<  ico.eo 

Ir( ICQ.EO 
I r  (  I C  9  •  E  Q 

eg  T?  107 

Ev>*CO  jE  (  °6 
CALL  XEXT 
JF(  IEP.NE 
IF(  ICO.EQ 
G9  T0  109 
ENCODE ( 96 
CALL  TEXT9(  1/ 
IF ( ICO  »E0 

go  T9  in 

ENCODE ( 96 
CAuL  TEXT 
IF<  IC9.E3 

Ir( ICO.EO 
IF(  IC^.EO 
[F( IC3.EQ 
IF(  IC9.EQ 

IF(  ICC.E0 


I  )  -V*SIV.(PSI  ) 
I )+V*CSS( DSI ) 

oT 

T*100.  ) 

T*100.  ) 

T*100*  ) 

0.  ) 

0«  ) 

H  1)102#105 


Uf 


TTOiTO^T3EA-^EA0j  ATVECOiCCC/TvEC 


< 


G6  T^  112 
E^C9DE<  q6 
CALL  TEXT 
IF( tCQ*EO 
39  TO  114 
ENC90EO6 


3)00 
15)  0 
22)0 
2?  )  3 
36)0 

5 


>30  TO 

TO  653  - 
T0  653 
1  TO  653 


101 


TO 
TO 
TO 


653 
653 
653 


(1, 


1/ 

I 

O)oo 

5)30 


ILAP  )  YC  /  X3,  ZZ/  DS  I  i  PH  I  /  T'-iET  A 


LAP/ 24/ 6/ 

T~'JT<  101) 

T6  108 


,3, IEP) 


ILK/ 


P°ESENiT 


10)  Oc 


ILAP  )  h2_j  Hi  /  H3/  44/  ~o*  ^6 
LAP/ 24, 8, 1/1^3/ IEP) 

°>    TO  110 


526/ 
(1*1 

2)33 
°  )  1a 
16)0 


3  0  )  0 

37)3 


ILAP 5 -31/ J30,-3  2/-33,^3^, h35 
LAP/ 24,9/ 1/1/3/ IEP) 

TO  601 

TO  601 
0  TO  601 

8  TO  601 

9  TO  601 
0  TO  601 


507 j  ILAP)L/Vz  ■• 
(1/ I  LAP/ 24, 14 j 

15)00  TO  113 

527/ ILAP)h7/HS 


>   H 


3/  111) 


H 9  »  -i  1 


-111 


CALL  TEXT8( 1/ lLAP/2^/16/ 1/ 1/3/ IER) 
!*"(  ICO.EQ.  20)00  T?  H5 
GO  T0  116 
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DATE 


115 

116 


6C4 
117 

113 

119 

120 
1  21 


602 


923 


ENC9D 

CALL 

if(  rc 
if<  t: 

IF(  IC 

Ir(  IC 

if(  rc 

IF(  IC 

ce  Te 

ENC9C 
CALL 
IF(  IC 


E( 

TC 


e  t 
l 

E( 

TE 


1 
EviC°DE{ 

TE 
•ft . 


CALL 

if<  rc 

G^  "ra 

EMC^D 

CALL 
I F  (  I C 
Ir  (  IC 
IF(  IC 
Ir  (  IC 
Ir(IC 
Ir<  Tc 
G9  T9 
EMC5P 
•n 


i 

i 

E( 

TC 

9. 
8. 
5. 

e. 


96,528,  I  la3  JH36/H37  ,H38,H39jH40  >-^t 
XT8(1/ ILAP*24* 17* 1*1*3* IER) 
EQ»3)G9    T8    604 
E0«10)G8    T3    604 
EQtl7)G8    T9    fcO^ 
EQ. 24)08    T8    604 
EQ.31)G8    T3    604 
•EG.38)GS    T3    604 
17 

96/ 511  j  ILAP)UD8T*V09T/'wD9T*RD9T/3D9T/QD8T 
XT8(  1/  ILAP,2<*,22>  1*  1/3/  IE??) 
♦25)09    T9    US 

,529,  ILAP)Hl3iMi4**l5j416/H17*  ms 
9( 1* ILAP/24/24^1/1j3j IE9) 

•  3  0 ) G 8    T  9     120 

6,53C/lLA=)i-!i*2/->43/-<+4,-!45,Hi<.6/-^7 
8(  1/  ILA=,2^,25,  1/  1,3,  IE9J 
.4)G5  "9  602 
T9  602 
602 
602 


EQ 

19 
96 
XT 
EQ 
21 
96 

X 


.11)09 
•13)09 
.25)09 
.32)09 
.39) 09 


T9 
T9 

TS 
T5 


602 

602 


9 
E( 


23 

96,516,  lLAP)T,Z>3ELT^IME,A=>ITc'-i,^EA0,  arATE,3ETA,a^9_L*S  =  EE 


CALL  TEXT9< l/lLAP/24,30* 1*1*3/ IE3) 

IF ( A101.GT.-C.90 ) 39  Ta  92^ 
u  =  0  • 
v  =  0 » 


RsO, 
P*Q, 

0  =  0. 
X  0  ■  0  • 

v  C  =  0  • 

zz=o. 

=  SI=C 
PHI=0 
THETA 
LD:?T  = 
M  0  8  "•  = 
n09t  = 

*oeT  = 

P09^  = 


TIME* 
H 1  ■  0  • 
H2  =  0- 

h3  =  0. 
h4sQt 


=  C 
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DATE 


H5  =  n. 
H6-C 
H7»0. 
H  3  *  0  • 

H9*C  • 
H10  =  C 
H11«G 
H12  =  0 
H13  =  0 
H 1  4  a  C 
H15  =  0 
H 1 6-»  0 
H 1  7  =  0 
H18  =  0 
H30*Q 
H31*0 
H32  =  0 
h33=0 
h34  =  C 
h35  =  C 
H36  =  C 
h37  =  0 
H33  =  0 
m39  =  0 
M40  =  Q 

h^a  =  D 

H  4  3  =  0 
H  **  4  s  G 

H4b  =  0 

H46  =  C 

^7  =  : 
rftNT 


c  \  «« 


60  3 


TO 


J24  CONTINUE 

IF(  IC5.EQ.5)_ 
IF( rce.EQ.12)G8 
lr(  !Ce.EQ«13)G9 

IF<  Ice«EQ. 26)39  TS  603 


03 
'3  603 


IF(  lce*EQ.33)G9  T?  603 
IF(  ICS.EG-^Q] 36  15  603 
35  TG  610 

603  C9NTINUE 

5  3   CONTINUE 

YiYs«45»SIN (PHI  ) 
XlXA=(2.5*R)+XlX 
Xlls2«5*R 

V1Y  a  =  Tw,ETA-y  1  v 
X1X3=(2»5*R)*X1X 
Yl  YE'  =  THETA  +  vi  y 

*(0« 1*3  IN (PHJ )  ) 


XlXC*(2«5*R)+( 0*1*3 IN(PHi ) ) 
Y1YCSTHETA+(C1*CSS(PMP  ) 
ICdA=ICSA+l-H 1 .♦U)+(10»*( IC5A/514- ) ) 
Yl  Y0= ( IC6A/51*. ) *( -1 .0)  +  ( ThETA ) -0.1  + 
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-q 


E5,n 


51 


55 


YlYE  = 
xlx^  = 
XlXE  = 
XlX^  = 
X 1  X  G  = 
Y1YP  = 
Y1YL  = 
Ha  (Ot 
ril=H* 
YlYEE 
REARA 
REARA 
REARS 
REA??6 
REAPC 
REAPC 
^EA?^ 

rear:, 
sesx  = 

SESY  = 
[F(Y1 

icsa  = 
c?nti 

I V  I E a 

I V  I  E  a 
IV  I  E  a 
IVIEa 

ivie* 

IVlE-v 
I V  I  E  a 
I  V  I  E  .v 
I V  I  E  M 
IF(R. 
IF(R. 
IV  IE* 
I V  I E  * 
IV  IE.. 
I V  I  E  a 
IV  IE.-. 
3?  T8 
LVIEW 
[VIEW 
IV  IE* 
IV  IE* 
IV  IE* 
38  Tg 
I V  I  E  A 
IV  IE* 
[VIEW 
I V  I  E  ri 
I  V  I  E  n 
38  T8 
C8NTI 


f  I 
(X 
XI 
VI 
(X 
vi 
vl 
2/ 
SI 

=  Y 

1  = 

2  = 

1  = 

2  = 

1  = 

2  = 
1  = 
c  = 
0. 
^  • 

XE 

n 

MJ! 

A 
A 
A 
A 
A 
A 
A 
A 
A 

LT 

A 
A 
A 
A 
A 

A 
A 
A 
A 
A 


C3A/ 

1X5- 
X5-C 
XA+C 
1  X  A  + 
YA-0 
YB-0 
C8S( 
\|  (  put 

1YE- 

-0.5 

-0.5 

-0.5 

-0-6 

C5 

-0.5 

0«5 

-0.6 

75+( 

75+( 


)  +  (THETA  )  -0.14 
A/514. )*0»2) 


514,  ) *( -1  .0 

0.05) -( I IC8 

.05 

•  05 

0*05) +< ( I C9 A/51  4. ) *0«2) 

.14 

.14 

PH  I  )  ) 

I  ) 

(  I  C8 A/2700. ) 

+  vl  Y 

4  +  Y  1Y 

-vly 

4-Y1Y 
YC/5000.  ) 

XC/5C0C  ) 
-0.6)38  T8  59 


2) 
3) 

5) 


)  = 

)  = 
)  = 

.0 


,0) 


IUEA 
IPAC 
IPaC 
IDAC 
IPAC 
IPAC 

IPAC 
I?AC 

IPAC 
2)38 
02)3 
■  IPA 


)*IPA 


2) 

3) 

4  ) 


=  I  =  A 
slPA 
■  IPA 


0(3* 
KCRE 
<(X1 
<(X1 

<(RE 
K(RE 
K(RE 
X(RE 
'<(RE 

T8 
8  T8 
C*(X 
C<(X 
C'<(X 
C<(X 
C'<  ( o 


10) 

ARA1 
X2,v 

XAj-Y 
ARC1 

A~A1 
A  RBI 
ARD1 

A^Cl 

51 
1XD* 
IXEj 
UFj 

Ixs, 
•  c*o 


>R£ARA2/0) 

lv3,  1  ) 
1YA#  1  ) 
,REARC2,  1) 
/REARA2,  l  ) 
,REAR32i 1 ) 
/REAR02i  1  ) 
>3EAPC2^  1) 


YIYEE^O) 
YiYE,  1  ) 
YlYEiC) 
YlYEEil  ) 


0 ) =I°ac< (XlXAj 

1)=IPAC<(X1XAj 

2)  =  I=3AC<(:?EAR3 

3)  =  ICAC< (XlxF* 
<0*IPAC<(X1xGj 


A  (  10)  »I°ac<  (X1X5# 
A(  11  )*lPAC'<(XlX3i 

A  (  12)  SIDAC'<  (^EaRB 

A(  13)  =  I°Ac<(XlxEi 
A(  14)  =IPAC'<  (XlX0# 

55 
MUE 


Y1YA/0) 

yiYRj  l ) 

1jREARD2j  1 

Y1YE>C! 
YlYEE; 1  ) 

Y1Y3.#0) 
rlYL;  1) 
1/REARB2;  1 

YlYEiC) 
Y 1 Y  E  E 1 1  ) 
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DATE 


610 


I V I Ew A { 15 ) » I P ACK  C  SESX j  SESY/ 0  J 

IV  IF  -,A(  16  )  =IPACK(SESX*SESYi  1  5 

IVIEWAt 17)»0 

CALL  GRAPHS (2, [V IEWA#17j2* IER ) 

ip ( rce.LE.^ojGa  t?  101 

IC8  =  0 

G9  re  lei 
end 
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